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ABSTRACT 

We present Spitzer imaging of the metal-deficient (Z ~30% Z ) Local Group dwarf galaxy NGC 6822. 
On spatial scales of ^130 pc, we study the nature of IR, Ha, H I, and radio continuum emission. 
Nebular emission strength correlates with IR surface brightness; however, roughly half of the IR 
emission is associated with diffuse regions not luminous at Ha (as found in previous studies). The 
global ratio of dust to H I gas in the ISM, while uncertain at the factor of ~2 level, is ~25 times lower 
than the global values derived for spiral galaxies using similar modeling techniques; localized ratios of 
dust to H I gas are about a factor of five higher than the global value in NGC 6822. There are strong 
variations (factors of ~ 10) in the relative ratios of Ha and IR flux throughout the central disk; the 
low dust content of NGC 6822 is likely responsible for the different Ha /IR ratios compared to those 
found in more metal-rich environments. The Ha and IR emission is associated with high-column 
density (i>10 21 cm" 2 ) neutral gas. Increases in IR surface brightness appear to be affected by both 
increased radiation field strength and increased local gas density. Individual regions and the galaxy as 
a whole fall within the observed scatter of recent high-resolution studies of the radio-far IR correlation 
in nearby spiral galaxies; this is likely the result of depleted radio and far-IR emission strengths in 
the ISM of this dwarf galaxy. 

Subject headings: galaxies: dwarf — galaxies: irregular — galaxies: ISM — galaxies: individual 
(NGC 6822) — infrared: galaxies 



1. INTRODUCTION 

Nearby dwarf galaxies provide a unique opportunity to 
resolve the interaction of the multiple phases of the in- 
terstellar medium (ISM) . Star formation in these systems 
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is more susceptible to small-scale processes (e.g., turbu- 
lence, feedback) than to the galaxy-wide perturbations 
(i.e., spiral density waves) that are common in more mas- 
sive galaxies. With their typically low nebular metallic- 
ities, star-forming dwarf galaxies serve as fiducial exam- 
ples for comparison with the high-redshift star-forming 
systems predicted i n the cold dark matter paradig m (e.g., 
Babul & FergusonllTflM lE7j^lT997t IMa,tedlT99rf) . 

The detailed processes of heating and cooling of the 
ISM are fundamental to our understanding of the na- 
ture of galaxy evolution. An environmentally-dependent 
fraction of UV photons produced by massive stars will 
be absorbed by dust and gas and re-radiated in the far- 
infrared (FIR) or in nebular emission lines. This fraction 
will depend on various factors, including the local dust- 
to-gas ratio (related to the metallicity and hence the dust 
content) and the porosity of the local ISM. In cases where 
the dust and stars are well mixed, and where the extinc- 
tion is moderate, there will exist a correlation between 
the fluxes in various wavelength bands, including the ul- 
traviolet (UV), Ha, FIR, and radio continuum. Nearby 
dwarf galaxies offer a unique opportunity to study one of 
the regimes where these correlations may begin to break 
down. 

NGC 6822 is a Local Group (D = 490 ± 40 kpc; Mateo 
1998) dwarf irr egular galaxy w ith a nebular abundance 
Z ~ 30% Z Q ijSkillman et alJ fl989t iLee et alJ l2006al). 
The total optical luminosity is ~9.4 x 10 7 L i|Mateol 
1998), while the total H I mass is 1.34 x 10 8 M Q (see 
detailed discussion in § 12. 21 and de Blok & Walter 2006). 
The system contains distributed low-level star forma- 
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Fig. 1. — DSS image of NGC 6822, overlaid with a contour of H I coulmn density at the 2 X 10 20 cm -2 level. The H I disk is much larger 
than t he high surface brightness optical body, which was the target of our IRAC and MIPS observations. However, de Blok & Walter 
(2006) find an extended stellar population that exceeds the size of the H I disk (see further discussion in £j HI. 



tion throughout the disk as well as actively star-forming 
regions wit h associ a ted H a emission which were first 
studied by Hubble (1925). The star formation has 
been quiescent over the past few Gyr, without ma- 
jor p e aks or dips in the absolute rate ijWvdeJ 120011. 
2003). Ide Blok fc Walterl l|2006j) derive a total Ha flux 
of 2 x 10 39 ergs -1 , corresponding to a global Ha-based 
star formatio n rate (SFR) of ~ 0.015 M Q yr 1 using the 
calibration of lKennicuttJ l|1998|) . 

NGC 6822 boasts an unusual H I distribution; in the 
extended H I disk (see Figure lies one of the largest 
known ho les in the ISM of a dwa rf galaxy, first dis- 
cussed bv Ide Blok fc Waited (|2000^ . Ide Blok fc Walterl 
( 2006) argue that this hole is likely caused by the com- 
bined effects of stellar evolution. That study also re- 
veals the presence of a significant blue star population 
in the outer H I disk, and a lower surface brightness, 
spheroidal red stellar population extending beyond the 
H I distribution. The extended H I disk may include a 
dwarf companion galaxy in the northwest region (note 
that interactions with this system may have played a 
part in triggering some of the recent star formation in the 
galaxy; see de Blok & Walter 2000). The well-sampled 
rotation curve sho ws that the system is highly dark- 
matter dominated ijWeldrake et all 12003). The proxim- 



ity of NGC 6822 means that it is one of only a few places 
where the small-scale (~100 pc) structure of the ISM can 
be studied at multiple wavelengths. 

Previous infrared studies of NGC 6822 have exam- 
ined the heating and cooling mechanisms in the ISM. 
Using InfraRed Astr onomical Satellite (IRAS) data, 
Gall agher et al.l l)1991f) found a variable Ha/FIR ratio 
within the disk, with roughly 50% of the FIR emission 
arising from regions that are luminous at Ha. The re- 
maining fraction of diffuse IR flux was attributed to a 
low optical depth (i.e., low dust content) in the ISM. 
Thus, regions of strong IR emission correspond to com- 
paratively dusty regions that also contain a y ounger stel- 
lar po pulation. A second IRAS study by llsrael et all 
( 1996) found that three of the bright IR peaks were as- 
sociated with luminous H II regions, while two equally 
bright sources were unassociated with obvious sources of 
ionizing radiation. These authors confir m the prominent 
diffuse dust continuum emission found bvfGallaghc r et all 
(|1991D . Further, they find, using single-temperature 
blackbody fits to the far-IR data, a dust-to-gas ratio 
(1.4 x 10~ 4 ) that is much lower than typical values seen 
in more metal-rich systems such as the Milky Way (av- 
erage dust-to-gas ratios of ~0. 006-0. 01; Sodroski etal. 
[1997], Li [2004]). NGC 6822 was observed with the In- 
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Fig. 2.— H I (a) and IRAC (b) images of NGC 6822. (a) shows the total H I column density distribution (de Blok & Walter 2000), while 
(b) shows the underlying stellar and warm dust components. The three-color image shows the 3.6 fim band as blue, the 4.5 fira band as 
green, and the 8 /im band as red; regions of hot dust emission, indicative of active star formation, appear as diffuse regions of red emission. 
The box in (a) shows the approximate field of view shown in (b). 



frared Satellite Observatory at 6.75 and 15 /m i, but only 
the most luminous H II re gions were detected ( Dal e et all 
2000| iHunter et afll200H) . 

In this paper, we present Spitzer Infrared Array 
Camera (IRAC) and Multiband I maging Photometer 
for Syitzer (MIPS) imaging (see I Werner et al.1 120041: 
iFazio et all 120041: iRieke et all l2004j) of the high optical 
surface brightness component (i.e., the inner portion of 
the H I disk) of NGC 6822 (see Figure for a compar- 
ison of the field of view imaged and the size of the H I 
disk). These data were obtained as part of the Spitzer 
Infrared Nearby Galaxies Survey (SINGS; see Kennicutt 
etal. 2003). As the nearest galaxy in the SINGS sample, 
these data probe the nature of mid- and far-IR emission 
at unprecedented sensitivity and resolution (~130 pc; see 
further discussion below). We compare these Spitzer data 
to sensitive Ha, H I and radio continuum imaging in or- 
der to study the relation between H I surface density 
and the locations and intensities of star formation and 
IR emission. 

2. OBSERVATIONS AND DATA REDUCTION 

2.1. Spitzer Imaging 

For an overview o f the S INGS observational strategies, 
see lKennicutt etal.1 l)2003j) . NGC 6822 was observed for 
122 minutes in IRAC mosaicing mode on 2004, October 
7 and 11; the separate visits to the target ensure proper 
removal of image artifacts and asteroids. The observa- 
tions cover the optical extent of the galaxy to the 
level (see Figures □ and EJ). Ide Blok fc Walteil (2006) 
find that the stellar population extends throughout the 
entire H I distribution; our IRAC data remain insensitive 
to stars in the outer H I disk of NGC 6822. The SINGS 
IRAC pipeline processes the basic calibrated data (BCD) 
images, producing mosaics with pixel scales of 0.75" and 
PSF FWHM values of 1.66", 1.72", 1 88", and 1 . 98" a t 
3.6, 4.5, 5.8 and 8 /jm, respectively l)Fazio et al.l 12004): 
note that 1" = 2.4 pc at the adopted distance of 490 kpc. 
Flux levels are uncertain at the ~ 10% level, primarily 



due to systematic effects in the calibration process. 

MIPS scan mapping mode observations were obtained 
on 2004, September 24 and 25, for a total of 162.1 
minut es. The MIP S Instrument Team Data Analysis 
Tool l)Gordon et al.ll2005l) was used to process the BCD 
files, producing mosaic images with pixel scales of 0.75", 
3.0", and 6.0" at 24, 70 and 160 /im, respectively. Sys- 
tematic uncertainties (e.g., detector nonlinearities, time- 
dependent responsivity variations, background removal, 
etc.) limit the absolute flux calibration to ~ 10% in the 
24 /im band and to ~ 20% in the 70 and 160 /im bands 
(note that the absolute calibration changes slightly de- 
pending on which pipeline was used to produce the BCD- 
level data 14 ; the S10 calibration is applied to these data). 
The FWHMs of the PSFs are 6", 18", and 40" at 24, 70, 
and 160 /im, respectively. All flux densities were ex- 
tracted from the broadband images after convolution to 
the 160 /xm beam. We use convolution kernels that con- 
vert an input PSF into a lower resolution output PSF 
using the ratio of Fourier transforms of the output to 
input PSFs. High frequency noise in the input PSF is 
suppressed when these kernels are created. See K. D. 
Gordon et al. (2006, in preparation) for details. 

At the Galactic latitude of NGC 6822 (-14.8°), the 
foreground stellar contamination is high. The 3.6/8 /im 
and 8/24 /im colors were examined in order to remove 
the brightest foreground sources, which will have the 
most pronounced effect on derived flux densities. Less 
luminous sources require a more detailed treatment; a 
first-order correction can be estimated by computing 
the strength of foreground star emission in regions that 
are clearly unassociated with the stellar component of 
NGC 6822, and sc aling these to th e sizes of the regions 
under study (e.s. ALee et al.ll2006b|) . However, the prox- 
imity of NGC 6822 makes it large on the sky; the IRAC 
maps do not include a large overscan region, making this 

14 See the MIPS data analysis handbook for details; 
http: / / ssc. spitzer . caltech.edu /mips 
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Fig. 3. — 2MASS J (a), 3.6 and 8 (im (b, c), 24, 70, and 160 fim (d, e, f), Ha (g), 20cm radio continuum (h), and H I spectral line 
(i) images of NGC6822. All image scales are linear, except for the Ha image (g), which is shown in a logarithmic stretch to highlight 
the low-surface brightness ionized gas that permeates much of the central disk of the system. Beam sizes are shown as boxed ellipses for 
the MIPS and radio data, in the bottom right corners. Overlaid are the sizes and locations of the smaller apertures [diameter = 114" ~ 
3 X FWHM (160 fim) ~ 270 pc] discussed in the text; these regions are labeled in panel (a). Note that the regions that encompass the 
major IR emission complexes (regions 3, 7, 8, 11, and 14; see also Table^ have also been analyzed using apertures with four times the area. 
The 8 fim image does not cover the full field of view shown in panel (c); since the 4.5 /an band has the same orientation, flux densities are 
not available for region 13 in either of these bandpasses. Note that the radio continuum image in panel (h) shows only the apertures where 
20 cm emission is detected; numerous background sources exist within this field of view (see further discussion in § 13.31 , The aperture sizes 
are large enough to cover the bulk of the extended long-wavelcngth PSFs, minimizing aperture correction effects (see further discussion in 
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techiiique susceptible to small-number statistics. We dis- 
cuss the effects of foreground contamination in more de- 
tail in § 13 Background sources are more difficult to 
identify (though a first-order estimate of their number 
density is available from some of the plots shown in this 
work; see further discussion in § 13. 2. 1(1 . and no correction 
for these sources has been explicitly applied. 



Foreground Milky Way cirrus conta mination is non- 
neglig ible in this region of the sky (e.g., iGallaeher et alJ 
1991); given the extent of the galaxy, there will be varia- 
tions in the cirrus contamination levels within the system 
that we are unable to remove. In an attempt to subtract 
the smooth diffuse component, high-surface brightness 
emission that is clearly associated with the galaxy (by 
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TABLE 1 
Infrared Sources in NGC 6822 



Identification 


a (J2000) 


6 (J2000) 


Previous I.D. 


References 


NGC 6822- 


1 


19-44-31 28 


— 14-47-4R fiQ 


196 1 


3 


NGC 6822- 


2 




— J.*±. 'X f J. ■ 1U. IU 


Huhhlp TI 

1 1 LI U Ult_ _L± 


I 


NGC 6822- 


3 


1 QfS 


— 14-49-10 71 




1 9 a 3 


NGC 6822- 


i 


1 9-44-37 89 


— 14-51 -37 7Q 


196 3 


3 


NGC 6822- 


5 


1 fiQ 
_i_ y . i i i ±. tj.uc' 


— 14-^fl-1Q £fi 


196 3 


3 


NGC 6822- 


6 


19:44:46.18 


-14:46:19.89 






NGC 6822- 


7 


19:44:49.47 


-14:52:40.93 


Hubble IV, G91 1, 196 4 


1, 2, 3 


NGC 6822- 


-8 


19:44:52.81 


-14:43:04.97 


Hubble V, G91 3, 196 5 


1, 2, 3 


NGC 6822- 


-9 


19:44:55.69 


-14:50:17.00 






NGC 6822- 


10 


19:44:57.36 


-14:44:50.01 






NGC 6822- 


11 


19:44:58.59 


-14:47:47.02 


Hubble VI/VII, G91 2, 196 6 


1, 2, 3 


NGC 6822- 


12 


19:44:59.82 


-14:54:41.03 




NGC 6822- 


13 


19:45:04.99 


-14:57:29.07 


196 7 


3 


NGC 6822- 


11 


19:45:05.01 


-14:43:17.07 


Hubble X, G91 3, 196 8 


1, 2, 3 


NGC 6822- 


15 


19:45:09.35 


-14:48:44.10 


Hubble IX, G91 2, 196 9 


1, 2, 3 


NGC 6822- 


16 


19:45:16.68 


-14:45:24.26 







References. — 1 - IHubbld fi925ft : 2 - IGallagher"etaLl IT99ll) : 3 - Hsrael etafl Il996ft 

a The sources in Gallagher ctal. (1991) are not defined by low-resolution (1.5' X 5') IRAS 
imaging; the cross-identifications here correspond to only the highest-surface brightness 
IR peaks in the IR or Ha images. 



comparison with images at multiple wavelengths) was 
masked, leaving an estimate of the (spatially variable) 
background cirrus emission. This background informa- 
tion was smoothed by a 6' boxcar function and then in- 
terpolated over the area of the galaxy. This smoothed 
cirrus component was then subtracted from the images 
at 24, 70 and 160 /im. Comparison of flux measure- 
ments before and after the cirrus removal shows that 
this correction is of order ~10%. The subtracted com- 
ponents at 70 and 160 are in rough agreement (in 
terms of both morphology and surface brightness) with 
the lower-reso l ution COBE/DIRBE maps presented by 
iSchleeel et alJ l)1998ft . Given the systematic uncertain- 
ties in the long-wavelength bands (see above) and the 
potential contamination by foreground cirrus emission, 
and assuming a random error component of up to ~ 10%, 
these corrections imply a total error budget of ~ 15% at 
24 fxm, and ~ 25% at 70 and 160 fim. 

2.2. H I Spectral Line and Radio Continuum Imaging 

The H I moment zero (representing integrated col- 
umn density) image pre sented here, first published in 
Ide Blok fc Waltel (2000), was obtained with the Aus- 
tralia Telescope Compact Array 15 (ATCA) in 1999 and 
2000, using the 375, 750D, 1.5A, 6A and 6D configura- 
tions. Fifteen 12-hour synthesis observations were taken 
in mosaicing mode, with the spectral line correlator pro- 
viding 0.8 kms -1 channel separation. 

To assure accurate total column density calibration, 
the intcrferometric data were combined with single- 
dish observations taken with the 64-m Parkes telescope. 
These data use the same channel separation as the inter- 
fcromctric observations. The "zero-spacing" correction 
was applied in the MIRIAD 16 environment. Conditional 

15 The Australia Telescope is funded by the Commonwealth of 
Australia for operation as a National Facility managed by the Com- 
monwealth Scientific and Industrial Research Organisation. 

16 See http://www.atnf.csiro.au/computing/software/miriad 



transfer functions (i.e., blanking) were applied to the fi- 
nal cube to create the integrated column density image 
presented here. The image has a beam size of 42.4" x 12" 
(1 JyBm" 1 = 1190.5 K). The final cube that was im- 
aged to produce the moment zero map has a 5 a column 
density sensitivity of 1.6 x 10 19 cm -2 ; the peak H I col- 
umn density is ~ 3 x 10 21 cm~ 2 , and the total H I mass 
is 1.34 x 10 s M Q at the adopted distance. More detailed 
discussions of the H I data, analy sis an d interpretation 
can be foun d inlde Blok fc Walter] (12000ft , iWeldrake etall 
( 2003), and de Blo k fc WalteT i|2b06'). 

During the acquisition of the ATCA H I data, a second 
frequency band was simultaneously observed; centered 
at 1.38 GHz, this 128 MHz bandwidth channel provides 
a sensitive L-band radio continuum map of NGC 6822. 
Note that these data have no zero-spacing correction ap- 
plied; the presence of comparatively strong background 
sources and the weak nature of the intrinsic radio con- 
tinuum in NGC 6822 make integrated total flux densities 
difficult to determine. Single-dish observations at other 
frequencies l|Klein et alll983|:lKlein fc Graevell986ft sug- 
gest an integrated L-band flux density of 100 ± 25 mJy 
ijlsrael et al.lll996j) . Comparing with our interferometric 
flux density measurements (see § below) suggests that 
the high surface brightness radio components account for 
^70% of the (uncertain) total L-band flux density. The 
remainder may exist as a diffuse synchrotron component 
(see more detailed discussion on the thermal/nonthermal 
decomposition within the galaxy in § I3.3fl . The radio 
continuum data use the same basic reductions as for the 
spectral line data, though the final image has a smaller 
beam size of 30" x 8"; the rms noise is ~ 90 /iJyBm -1 . 
Note that the 20 cm image presented in Figure|3Jh) shows 
only the apertures with radio detections. Further, there 
are numerous background radio sources in the field (see 
also the discussion in § 13. 3|) . 

2.3. Deep Ha Imaging 
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TABLE 2 

Observed Quantities in NGC 6822 a 



T-> ■ b 

Region 


Ha 

(ergs -1 cm -2 ) 


3.6 (im c 
(mJy) 


4.5 /im 1 
(mJy) 


5.8 /im c 

(mJy) 


8.0 fj,m c 

(mJy) 


24 fim 
(mJy) 


70 /im 
(mJy) 


160 
(mJy) 


20 cm 
(mJy) 


H I 

(Jykms -1 ) 








Aperture 


radii — o t - 


= 3 x (160 Atm PSF FWHM ^ 130 pc) 








1 


(9.6 ± l.» )hj-Lo 


24 ± 3 


16 ± 2 


16 ± 2 


28±4 


100 ±16 


560 ±140 


o M nn _L con 
zlOO ± ooO 


^1.0 ± 0.3 


12.9 ± 1.3 


z 


( OK -1- C.M7 1 Q 

(zo ± )iii-Lo 


1 7 4- O 
1 i ± Z 


1 O 4- O 
1Z It z 


p. 4- n « 
it U.O 


15 ±2 


20±3 


470 ± 120 


01 nn 4- ^Qn 
zlUU it OoU 


t Q 4- fl A 


1 O O 4- 1 O 
1Z.Z It l.Z 


3 


1 1 i n _l_ on MP i q 
(110 ± z0)r.-13 


17 ± 2 


14 ± 2 


11 ± 2 


27 ±4 


110 ±17 


1700 ±430 


'J7AO _L (19(1 

o /00 ± 930 


12 ± 1.2 


12.0 ± 1.2 


A 


(z4 it Jrij-lo 


11 4- A 


01 4- 1 
Zl It o 


on 4- Q 
ZU It o 


26 ±3 


32 ±5 


630 ±160 


01 nn 4- ^Qn 
zlUU it OoU 


Q 1 4- n A 
o. 1 It U.4 


1 a a 4- 1 k 

14. o ± l.O 


5 


(5.0±1)E-13 


43 ±5 


29 ±3 


26 ±3 


27 ±3 


81 ±13 


600 ± 150 


1600 ± 400 


<;i.0±0.3 


11.0 ± 1.1 


6 


(14±3)E-13 


55 ±6 


38 ±4 


34 ±4 


33±4 


17 ±3 


520 ± 130 


2000 ± 500 


^1.0±0.3 


10.1 ± 1.0 


7 


(24±5)E-13 


42 ±5 


29 ±3 


49 ±6 


84 ± 10 


260 ±40 


2300 ± 580 


4900 ± 1300 


6.4 ±0.7 


21.8±2.2 


8 


(120±25)E-13 


49 ±6 


39 ±4 


65 ±8 


130 ± 14 


780 ± 120 


4300 ±1100 


6800 ± 1700 


21 ±2.1 


14.0 ±1.4 


9 


(4.6±0.9)E-13 


74 ±8 


49 ±4 


47 ±5 


41 ±4 


23±4 


900 ± 230 


2100 ±530 


£1.0±0.3 


13.7±1.4 


10 


(4.8 ± 1.0)E-13 


62 ±7 


43 ±4 


44 ±5 


41 ±5 


31 ±5 


900 ± 230 


1700 ± 430 


^1.0±0.3 


9.1±0.9 


11 


(13±3)E-13 


92 ±10 


60 ±7 


76 ±9 


100 ± 12 
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The Ha image presented here was obtained with the 
2.5m Isaac Newton Telescope, using three pointings 
with the Wide Field Camera to cover the entire H I 
disk of NGC 6822. The integration time was 80 min- 
utes per pointing; wide-filter R-band imaging was used 
to remove the continuum. Flux calibration was ob- 
tained by comparison with published fluxes fo r high- 
surface brightness regions in iHodge et all lll"98l []~989), 
after correction fo r foreg round reddening (E(B— V) = 
0.24; IGallart et al.lll996a|i . The limiting flux level is ~ 
3.7 x 10~ 18 ergs -1 cm~ 2 , roughly a factor of ten deeper 
than the images presented by Hodge and collaborators. 
For a detailed discussion of the Ha image handling, see 
Ide Blok fc Walterl ll200l . 

3. MULTIWAVELENGTH EMISSION IN NGC 6822 

Multiwavelength imaging of NGC 6822 reveals a wealth 
of detail and structure in the stellar and gaseous com- 
ponents. In Figure [3] we present images of the central 
region of the galaxy at nine dif ferent wavelengths. In 
the 2MASS J (Jarrctt ct al. 2003() and 3.6 /im bands [see 
Figures EJa) and^b)], the red stellar component of the 
galaxy is smoothly distributed, with few major stellar 
concentrations or clusters (compare to Wyder 2001, 2003 
and de Blok & Walter 2006 for more detailed views of the 
system at optical wavelengths) . Longward of 4 /xm, the 
spectral energy distribution (SED) of the stellar popula- 
tion has dropped off sufficiently that the warm dust emis- 
sion from active star formation regions dominates the 
galaxy morphology (see the 8 /xm image in Figure EJc). 



Moving into the FIR, Figures e, and f show that 
the active regions of star formation stimulate warm dust 
emission throughout the system. There is also very cool 
dust extending to the south of the major star formation 
regions (though still well within the optical and H I radii; 
see Figure OJ); comparing this emission with the Ha and 
H I images shown in Figures Elg) and |^i) leaves little 
doubt that it is associated w ith NGC 6822 and is n ot an 
image artifact. In deed, both iGallagher et alJ l)1991|) and 
llsrael et all <Tl 996 , recover this str ucture in IRAS maps. 
Further, Ide Blok & Walterl l)2003[) find that this region 
(spatially coincident with the rim of the H I hole) has 
the largest concentration of blue stars in NGC 6822. 

The Ha image presented in Figure Efg) shows that 
there is ionized gas throughout much of the central disk 
of the galaxy; previous observations and comparisons 
with IR imaging have concentrated on th e higher-surface 
brightness star format ion regions (e.g., Gal lagher et alJ 
119911 llsrael et al~lll996l, . Six of these regions correspond 
to detections in the ATCA radio continuum image shown 
in Figure|3Jh) . Note that there are numerous background 
sources in the radio continuum field of view; we have clas- 
sified radio "detections" as areas of > 5<r flux density, 
that have spatial counterparts in Ha and the 24, 70 and 
160 /im bands (note that we only show the apertures 
with radio detections in Figure |3fr) . Finally, the H I 
distribution (Figure |3) throughout the central region of 
the galaxy contains mostly high-surface brightness (i.e., 
H I columns J> 10 21 cm" 2 ) emission. There is clumping 
of the neutral gas surrounding various emission peaks in 
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Fig. 4.— The 24 (a, b, c), 70 (d, e, f) and 160 /jm images (g, h, i) of NGC 6822, with Ho (b, e, h) and H I (c, f, i) contours overlaid. The 
Ha contour is at the flux level of 3.8 X 10 — 18 ergs -1 cm -2 , highlighting low-level nebular emission throughout the disk; the H I contour is 
at the level of 10 21 cm -2 , highlighting the canonical surface density threshold for star formation. Note the striking correspondence between 
nebular, dust continuum, and H I emission. 



other wavebands. 

We study the panchromatic emission from 16 regions 
throughout the galaxy; these regions are shown as cir- 
cles overlaid on the various panels of Figure and their 
positions and cross-identifications are given in Table ^ 
The aperture locations were chosen primarily to probe 
the IR emission complexes. The aperture sizes and lo- 
cations represent a compromise between spatial resolu- 
tion, minimization of aperture effects, and diversity of 
regions that can be probed. The aperture diameter of 
114" (= 270 pc at the adopted distance) is ~ 3 x larger 
than the FWHM of the 160 fim PSF (and considerably 
larger than the FWHM of the 70 /xm, 24 iiin, and IRAC 



bands), which should alleviate serious aperture correc- 
tions to the measured flux densities (especially at 3.6, 
4.5, 24 and 70 /xm; extended-source aperture corrections 
at 5.8 and 8.0 are still uncertain, but may be as large as 
~30%). For example, an aperture-centered point source 
at 160 /Ltm would have an aperture correction of ~ 30%; 
this factor will be smaller for extended emission which 
fills the beam (as is the case for most of the 160 /jm 
sources in the galaxy; however, the size of the correction 
will depend on the IR surface brightness in the regions 
surrounding the apertures as well). This chosen aper- 
ture size is well-matched to most of the discrete dust 
emission regions in the galaxy, while still allowing us to 
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Fig. 5. — Surface brightness profiles of H I emission, overlaid with profiles of Ha emission (a), and 24, 70 and 160 [ira emission (b, c, d). 
The solid line and left-hand y-axis label in each panel correspond to the H I surface density in units of Mq pc — 2 (not inclination corrected); 
the dotted lines and right-hand y-axis labels correspond to the flux levels in Ha (in units of ergsec^ 1 cm" 2 ) or to IR surface brightness 
values (in units of MJysr" 1 ). These profiles were created using the same parameters as the rotation curve analysis of Wcldrakc ctal. 
(2003); the central location is (a, S) = 19:44:58.04, —14:49:18.9 (J2000; this location is between apertures 9 and 11 in Figure El, and an 
average position angle ~ 120° is used. The profiles integrate across the high-S/N (> 3<r) regions of the images shown in Figurcsl3l and |41 



resolve variations in the nature of the FIR emission by 
using a large number of apertures. Table|21shows the flux 
densities derived in each aperture from Ha to the radio 
regime. Note the footnotes to that table, which quan- 
tify the effects of foreground stellar contamination in the 
short-wavelength data. We also investigate the effects 
of expanding the apertures by a factor of 4 in area for 
the five brightest dust emission complexes (see further 
discussion in § 13. If) . Table |3] presents various quantities 
derived from the measurements. We derive an estimate 
of the 3- 1000 /xm total infrar ed flux (TIR) using equa- 
tion 4 of lDale fc Heloul l)2002|) . evaluated at the effective 
MIPS wavelengths: 

TIR = (10~ 14 )(19.74 • S* 24 + 3.23 • S 70 + 2.59 • S wo ) (1) 

where S24, S70 and S160 are the measured flux densities 



in units of Jy. Table also presents the Ha/TIR ratio, 
the dust-to-gas ratio (see more detailed discussion below, 
and also Tables 01 and 0J, and the ratio of TIR/radio 
luminosity as represented by the radio-FIR correlation 
coefficient "q" ; each of these quantities will be discussed 
in more detail below. 

3.1. Linking Heating Sources with the Dust Continuum 

3.1.1. Ha vs. IR Emission 

The major H 11 regions Hubble I /III, V, and X (re- 
gions 3, 8, and 14 in Figure respectively) are the 
strongest Ha sources in the galaxy. These regions are 
also luminous at 24, 70 and 160 /im. However, the rela- 
tive strengths of Ha and TIR emission vary considerably 
from one source to the next; the regions of strongest Ha 



9 



i — ;] — i — i — rr 




i i i i i i i i i islj i i i i i i i i i i 

1 10 100 1000 

Ha Flux (10~ 13 erg s _1 cm -2 ) 

Fig. 6. — Comparison of measured fluxes at Ha and 24 (an; black points indicate the 16 regions shown in Figure |3] red points show 
the five larger apertures used for the major IR peaks (see Table |2*| : the blue point shows photometry for the entire galaxy. Asymmetric 
error bars on points for regions 7 and 8 (shown in cyan) demonstrate the magnitude of the me asured extinctions at Ho, as inferred from 
comparing (thermal) radio continuum and Ha emission strengths (see discussion in § 13.1.11 . The 24 /im flux is calculated using the 
techniques discussed in Calzetti et al. (2005), where F24 = S24 x 1*24; values of S24 (in Jy) are taken from Table |5] and the frequency is 
evaluated at the wavelength of 24 fim. The relation derived by Calzetti etal. (2005) for optically thick H II regions in M51 is shown by the 
solid green line. The dotted lines show relations of direct proportionality, as labeled, for reference. 



and TIR emission are not always co-spatial. There are 5 
sources in the galaxy with TIR fluxes > 1.5 Jy (regions 
3, 7, 8, 11, and 14; see Table|3J); amongst these regions, 
the Ha/TIR ratio varies by a factor of ten, with region 
11 having strong TIR emission but weak Ha, and with 
region 3 having strong Ha compared to TIR. Expanding 
these apertures by a factor of 2 in radius shows the same 
trends, but with reduced extrema (a factor of 5 variation 
in the Ha/TIR ratio is found between regions 11 and 
3). Similar scatter is seen across the entire sample of 16 
regions (see Table 

We now examine what fraction of the observed vari- 
ations c an be attribu ted t o extinction effec ts. Recent 
work bv lWvderl lj2003|) and lLee etal] l|2006a|) shows that 
there is variable extinction throughout the disk, with 



some regions showing no optical extinction and others 
showing Ay values up to ~2 magnitudes. For strongly 
star- forming regions with low Ha/TIR ratios (e.g., re- 
gions 7 and 11), the Ha flux (and, correspondingly, the 
Ha/TIR ratio) may be up to 6 times larger than ob- 
served af ter correction for potentially strong internal ex- 
tinction ijLee et alJl2006a|) . It should be noted, however, 
that these values were extracted using 1.5"- wide spectro- 
scopic slits, and thus probe extinction on a much smaller 
spatial scale than we are sensitive to with our Spitzer 
data. Thus, while extinction can play a localized role in 
some regions, it cannot account entirely for the observed 
variation in the Ha/TIR ratio. 

A second measure of the extinction is available by com- 
paring Ha and radio continuum emission strengths. For 
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Fig. 7. — Histograms of H I column densities where high-surface brightness Ha and FIR emission are detected. The thick black lines 
correspond to flux levels of ~ 3.7 X 10 — 18 ergs -1 cm -2 (Ho) or to surface brightness levels of 4 MJy sr — 1 (IR), and are enumerated by the 
left-hand y-axis of each plot; the thin gray lines correspond to flux levels of ~ 1.4 X 10 — 17 ergs -1 cm -2 (Ha) or to surface brightness levels 
of 16 MJy sr -1 (IR), and are enumerated by the right-hand y-axis of each plot. The dotted line shows the histogram of H I column densities 
throughout the entire galaxy, normalized arbitrarily. These plots demonstrate that nebular and dust emission are strongly correlated with 
regions that are rich in H I gas. 



this comparison, we assume that the 20 cm flux densi- 
ties shown in Table [5] arise from thermal emission pro- 
cesses. This assumption is justified for most of the major 
H II regions in the galaxy, based on two lines of evi- 
dence. First, multi-frequ ency radio con tinuum data us- 
ing both int erfero metric |C ondorJ ll987f l and single-dish 
{Klein et al]ll983 iKlein fc Graevdll986l) observations as 
compiled in llsrael et all <|1996j) show that NGC6822 is 
extraordinarily weak in the nonthermal radio continuum 
- nearly all of the emission throughout the galaxy ap- 
pears to be of a thermal origin Secon d, applying the 
relations in lCaplan fc Deharvend 119861). an d assuming 
the T e values calculated in lLee et alJ l)2006aj) . we can use 
our observed (i.e., not extinction-corrected) Ha fluxes to 
predict lower limits for the strengths of thermal radio 
continuum emission. This exercise shows that nearly all 
regions have a small or (in most cases) negligible non- 
thermal component: only regions 7 (Hubble IV) and 8 
(Hubble V) show thermal fractions less than 90% (~50% 



and 70% in regions 7 and 8, respectively). This suggests 
that these two regions, both active in current star forma- 
tion (as evidenced by their strong Ha and IR emission), 
have produced sufficient numbers of SNe to accelerate 
the relativistic electrons that give rise to the apparent 
nonthermal components. 

We now calculate inferred reddening values, proceed- 
ing with the assumption that thermal emission domi- 
nates throughout the galaxy. Of the six regions detected 
in the radio continuum, only regions 7 and 8 show po- 
tentially strong inferred reddening values (Ah q <: 0.25 
magnitudes). Regions 2, 3, 4, and 14 are consistent with 
values of Ah q < 0.1 magnitude. Note that if regions 7 
and 8 do in fact contain a (relatively) strong synchrotron 
component, then the (smaller) thermal fraction of the de- 
tected flux densities puts a yet stronger constraint on the 
low inferred values of the reddening. The smaller inferred 
extinctions using this approach co mpared to the values 
derived from long-slit spectroscopy ijLee et alJl2006a() are 
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Fig. 8. — Pixel-by-pixel comparison of the strengths of H I and Ha (a), 24 /im (b), 70 /im (c), and 160 /an (d) emission. From these 
plots it is clear that both nebular and dust emission is found preferentially in regions of high H I column density. The 24 and 70 fim plots 
(b, c) are contaminated by foreground stars and background galaxies at low column densities; a first-order estimate of the severity of the 
contamination is available by comparing these panels. 



expected, since the latter data probe smaller-scale dust 
concentrations. 

To investigate the localized effects of photon absorp- 
tion and re-radiation in the IR, the central column of 
Figure 0] presents detailed comparisons of the Ha and 
IR morphologies. Here, a low-level Ha flux contour 
(3.7 x 10~ 18 ergs -1 cm -2 ) is overlaid on each of the 24, 
70 and 160 /jm images. It is immediately clear that the 
regions of high Ha and IR surface brightnesses are di- 
rectly linked, since the Ha contours surround the IR 
emission peaks. We are thus observing localized dust 
heating in the ISM (see also further discussion below). 
However, it is also interesting to note that a significant 
fraction of the integrated IR flux density of the galaxy is 
unassociated with regions that are luminous at Ha. Most 
Ha-bright regions are also luminous in the IR; one no- 
table exception is region 16 (see Figure yj, which shows 
(low-level) Ha emission in a clear shell morphology but 



which lacks strong associated dust emission. However, 
we treat region 16 with caution, since the local back- 
ground is not well-defined (see Table[5J and may be over- 
subtracted in the IR images. 

To discern the origin of the photons that give rise 
to IR emission throughout the galaxy, it is useful to 
quantify what fraction of the total IR luminosity of 
NGC 6822 arises from regions that are not luminous at 
Ha. The global flux densities in the long- wavelength 
bands, measured over the entire central disk region 
(S 24 = 2.51 ±0.50 Jy; S 70 = 53.2 ±15 Jy; S 160 = 
136.2 ±40 Jy; see Figure Eland Table appear to con- 
sist of discrete FIR peaks and a more extended "diffuse" 
component. Here we classify emission as diffuse if it is 
located outside the regions detected at Ha (shown as 
the 3.7 x 10~ 18 ergs -1 cm -2 contour in Figure^). The 
quantification of the fraction of the total FIR that arises 
in a diffuse component is dependent on the sensitivity 
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Log[S(70 / Lxm)/S(l60 /j.m)] Log[S(70 / am)/S(l60 /xm)] 

Fig. 9. — Ratios of IR to H I emission strength as a function of dust temperature: (a) shows that as the dust temperature increases (i.e., 
higher radiation field strengths), there is more 70 /an emission per unit H I mass; this suggests that local star formation is an important 
parameter in determining the characteristics of IR emission, (b) shows that cooler dust is less sensitive to heating; while a marginal trend 
exists for higher 160 fim fluxes per unit H I mass in higher-temperature environments, it is less evident than at 70 /im (as expected). 
Small-radius (57") apertures are shown as black points, while larger-radius apertures (114") are shown in red. Note that region 16 is not 
plotted, since the local background is not well-defined. 



of the Ha data to which we compare the IR images (or, 
equivalently, on the surface brightness below which we 
classify IR emission as "diffuse"). The spatial scales over 
which the diffuse emission are found are much larger than 
the PSFs at 160 and (especially) 70 /im (precluding PSF 
smearing as the only origin of the diffuse component). 
We find that > 60% of the integrated flux densities at 70 
and 160 /im arise from this diffuse component. While the 
24 /im emission is more strongly peaked in regions of ac- 
tive s tar formation fe.g.. lHelou et alJl2004t ICalzetti et alJ 
|2005|K the diffuse component still accounts for ~ 50% of 
the total 24 /im flux density. These results suggest that a 
substantial amount of the IR radiation may be powered 
by non-ionizing sources or may be the result of UV pho- 
tons escaping directly from star formation regions; this 
is discussed furth er in jj|H Simila r results were obtained 
for NG C 6822 bv IGallagherltafl JT991 and llsrael etaD 
( 1996); diffuse far-IR emission is also seen in other near by 
dwarf galaxies (e.g., NGC 55; Engelbracht etal. 2004). 

In FigurclSlwe present elliptical surface brightness pro- 
files of H I, Ha and IR emission. These profiles were 
derived us ing parameters from the H I rotation curve 
analys is of Ide Blok fc Waited IpOOOl ) and lWeldrake et all 
l|2003|) . Radially -averaged elliptical annuli with 12" 
thickness are integrated along the position angle of the 
galaxy's major axis [central position: a = 19:44:58.04, 
5 = -14:49:18.9 (J2000; this location is between aper- 
tures 9 and 11 in Figure [3J; position angle ~ 120°; see 
Weldrake et al. 2003 for details] . The resulting surface 
brightness profiles show that Ha, IR and H I peaks cor- 
relate remarkably well. In particular, the Ha and 24 /im 
profiles are very similar, tracing the local star formation 
rate (Helou etal. 2004; Calzetti etal. 2005). We discuss 
the surface brightness profiles in more detail in § 13.21 

3.1.2. Infrared vs. Optical Luminosities and Star 
Formation Rates 



The metal-poor ISM of NGC 6822 is an interesting en- 
vironment in which to compare the luminosities and im- 
plied SFRs from various methods. Using the global flux 
densities in the long-wav elength bands (see a bove) and 
the relations presented in lDale fc Heloul l|2002j) . the total 
IR flux of ~ 5.7 x 10~ 12 Wm~ 2 corresponds to L T ir ^ 
1.6 x 10 41 ergs" 1 at 490 kpc. The Ha-based SFR (sen- 
sitive to recent star formation over the last ~ 10 Myr) 
is ~ 0.016 M yr _1 , and the total recent SFR (over the 
last ~ 1 00-200 Myr) derived f rom the stellar population 
study of iGallart etal] (|1996b|) is -0.04 M Q yr" 1 . 

The 24 /im luminosity (L24), which has been shown 
to correlate well with regions of active star forma- 
tion in a range of envi ronments (e.g., iHelou et alJ 
12004 iCalzetti et all l2005|) , is an appropriate IR-based 
SFR indicator for dust-rich regions. We compare the 
I^/Ha ratio for NG C 6822 with what was observed by 
Cal zetti et all (120051) in th e more metal-rich spiral galaxy 
M51 l|Bresolin et allEffll . as follows: 

Log(L 2i ) = lMLog{L Ha ) - 0.06945 (2) 

where L24 and Ljj a are the observed luminosities at 24 
/im and Ha, respectively; the intrinsic Ha/Pa ratio is 
assumed to be 8.734. The central wavelength of 24.0 /im 
is used to convert the 24 /im flux density to a luminosity 
measurement [thus, L24 as defined above is not strictly a 
luminosity, but is rather equivalent to v ■ F„; we keep the 
above notation to allow direct comparison to the work of 
Calzetti et al. (2005)]. Using the total observed Ha lumi- 
nosity of the system (see Table [2J , this relation predicts 
a total 24 /im luminosity —3 times higher than the ob- 
served value. This confirms that metal-poor galaxies do 
not have enough metals to form dust in the same fraction 
as more metal-rich environments. 

To explore the relation between Ha and IR emission 
on a spatially resolved basis, we plot in Figure |H1 the 
Ha vs. 24 /im fluxes, and compare to the relation de- 
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TABLE 3 
Derived Properties of NGC 6822 



Region 


S(TIR) b 


F(Ha)/F(TIR) c H I Mass 


M Dust /M H i d 


qTiR 


. , c 

q7o 


q24 




(10- 13 Wm 2 ) 




(10 5 Mq) 


IP) 










Aperture radii = 57" 


~ 3 x (160 /urn PSF FWHM ~ 130 pc) 








I 


1.03 ±0.41 


9.6 ±4.3 


7.3 ±0.8 


0.0038 ±0.0021 








2 


0.73 ±0.29 


36 ± 16 


6.9 ±0.7 


0.0045 ±0.0024 


2.76 


2.15 


0.78 


3 


1.98 ±0.79 


55 ±24 


6.8 ±0.7 


0.0048 ±0.0026 


2.66 


2.15 


0.96 


1 


0.85 ±0.34 


28 ±13 


8.4 ±0.9 


0.0039 ±0.0021 


2.87 


2.31 


1.01 


5 


0.80 ±0.32 


6.3 ±2.8 


6.2 ±0.7 


0.0034 ±0.0018 








6 


0.76 ±0.30 


18 ±8 


5.7±0.7 


0.0046 ±0.0025 








7 


2.88 ± 1.15 


8.3 ± 3.8 


12.4± 1.2 


0.0037 ± 0.0020 


3.08 


2.56 


1.61 


8 


5.39 ±2.16 


23 ±10 


7.9 ±0.8 


0.0073 ±0.0039 


2.83 


2.31 


1.57 


9 


0.91 ±0.36 


5.1 ±2.2 


7.8 ±0.8 


0.0026 ±0.0014 








10 


0.69 ±0.28 


7.0±3.2 


5.2 ±0.6 


0.0020 ±0.0010 








11 


2.40 ±0.96 


5.4±2.5 


8.6 ±0.9 


0.0066 ±0.0035 








12 


0.52 ±0.21 


11 ±5 


6.2 ±0.7 


0.0022 ±0.0012 








13 


0.70 ±0.28 


9.7 ±4.4 


8.0 ±0.8 


0.0030 ±0.0016 








11 


2.42 ±0.97 


46 ±20 


9.2±1.0 


0.0025 ±0.0013 


2.75 


2.34 


1.24 


15 


1.07 ±0.43 


15 ±7 


6.0±0.6 


0.0035 ±0.0019 








16 f 


0.29 ±0.12 


59 ±28 


6.6±0.7 


0.00055 ±0.00029 










Aperture 


radii = 114" 


~6x (160 fan 


PSF FWHM) ~ 260 pc 








3 


4.0±1.6 


33 ±14 


23.3 ±2.4 


0.0050 ±0.0027 


2.97 


2.50 


1.19 


7 


5.0±2.0 


8.4±3.7 


36.4 ±3.7 


0.0030 ±0.0016 


3.01 


2.56 


1.40 


8 


7.8 ±3.1 


19 ±8.6 


25.6 ±2.6 


0.0040 ±0.0021 


2.99 


2.54 


1.62 


11 


5.7±2.3 


6.1 ±2.8 


28.4 ±2.9 


0.0048 ±0.0026 








14 


4.1 ±1.6 


32 ±15 


30.1 ±3.0 


0.0016 ±0.00087 


2.80 


2.46 


1.18 


Total Galaxy 


NGC 6822 


57 ±23 


12.3 ±5.5 


1340 ± 140 


0.00077 ± 0.00043 s 


3.34 


2.88 


1.55 



a Scc Tabic 1 for source coordinates and cross-identifications. 

^Calculated using the relations presented in ID ale' & Hclou 12002): uncertainties arc estimated at the ~ 40% 
level. 

c Calculated as the unitlcss ratio of Ha to TIR flux; ratios represent 1000 X this fraction, 
bee Tabic |H for dust mass deriv ations; the va lues shown here are calculated using dust masses derived from 
the models of fDraine fcTl I2001D . 1lT& Draind I2001ll2003) . and lDraine & Lil I2006fl . 

e c[TiR is the "q" parameter evaluated using an estimate of the total IR luminosity (log(S(TIR/Radio))); q24 
and q7o are the monochromatic "q" parameters, as studied in Murphy et al. (2006). 

Note that the local background surrounding region #16 in the MIPS images is very uncertain; see 
s This value uses the total H I mass from the entire system; using the H I mass contained in the matching 
aperture used for the global IR flux density extractions reduces the enclosed H I mass by ~ 10%. 



rived for M51. While there is an appreciable scatter, 
it is clear that regions within NGC 6822 follow a dif- 
ferent scaling between Ha and 24 /im luminosity than 
the regions in M51: the H II regions in NGC 6822 have 
much lower extinction than those in M51, as expected 
on the basis of the widely different metallicities (factor 
~10-15; see Bresolin etal. 2004). Note that extinction 
corrections at Ha will move points to the right in this 
plot, and cannot directly account for all of the observed 
dispersion (see further discussion in § I3.1.1|) ; we plot the 
largest detected reddening vector (0.75 mag, for aperture 
7 or Hubble IV) in Figure El to demonstrate the severity 
of (spatially variable) extinction at Ha. The scatter in 
this plot can then be interpreted as the effect of star 
fo rmation on th e loca l ISM . This expands the findings 
of iCannon et alJ ((20051 l2006ft . where the vigorously star- 
forming, metal-poor galaxies IC 2574 and NGC 1705 were 
also found to have different Ha/24 /im ratios than those 
found in more me tal-rich galaxies. Further, it reiterates 
the conclusions of iKennicuttj (TlJ98): the infrared lumi- 
nosity is a reliable SFR indicator only in metal-rich (i.e., 
dusty) environments, where a large fraction of the light 
from massive stars is reprocessed by dust to the infrared. 
We discuss these results further in § 0] 



3.2. Neutral Gas vs. Dust Emission 
3.2.1. Correlations Between H I and Dust 

The right column of Figure^Jshows the IR images over- 
laid with contours of H I column density at the 10 21 cm~ 2 
level. It is clear from this image and Table El that the 
central disk of NGC 6822 is rich in neutral gas and has 
relatively small variations in H I column density over the 
apertures in this study (factor of ~ 2 variations in flux). 
Comparing the locations of Ha, H I and dust emission in 
Figure^] it is also clear that nearly all locations of active 
star formation (as evidenced by high-surface brightness 
Ha emission) are associated with H I columns in excess 
of ~10 21 cm~ 2 . Note by comparison with Figures^and 
12 however, that the H I distribution is much more ex- 
tended than the high-surface brightness optical body and 
the areas with strong FIR emission. 

The surface brightness profiles of the FIR emission (see 
FigureEJ) show that the 24 /im emission is, in general, as- 
sociated with local H I column density maxima. Moving 
toward cooler dust emission, the 70 and 160 /im profiles 
show a much lower dynamic range than the 24 /xm pro- 
file; however, the peaks in the FIR nicely correlate with 
the locations of H I surface density maxima. For exam- 
ple, at radii of ~0.2, 1.0, and 1.6 kpc, the correlations 
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TABLE 4 
Dust Mass Derivations in NGC 6822 





• a 


rn b 

J- DUST 


Blackbody 




Dale et al. 


Li &: Drainc 






(k) 


Dust Mass b 




Dust Mass c 


Dust Mass d 






(10 3 M ) 




(10 3 M S ) 


(10 3 M ) 




Apcrt 


lug radii — 


57" = 3 x (160 ixra PSF FWHM £ 130 pc) 


I 




21 ± 2 


1.4±0.7 


2.03 


12.6±6.3 


2.81 ±1.40 


2 




21 i 2 


1.6 ±0.8 


2.88 


14.4 ±7.2 


3.11 ±1.55 


Q 





24 ± 2 


1.6 ±0.8 


2.31 


14.4 ±7.2 


3.27 ±1.63 


4 




22 i 2 


1.3 ±0.6 


2.75 


11.7 ±5.9 


3.26 ±1.63 


5 




23 ± 2 


0.82 ±0.41 


2.00 


7.4 ±3.7 


2.11 ±1.05 


6 




21 ±2 


1.4±0.7 


2.81 


12.6 ±6.3 


2.63 ±1.31 


7 




24 ±2 


2.1 ±1.1 


2.00 


18.9 ±9.5 


4.62 ±2.31 


8 




25 ±2 


2.3± 1.1 


1.72 


20.7 ±10.4 


5.79 ±2.89 


9 




23 ±2 


0.97 ±0.49 


2.62 


8.7 ±4.4 


2.04 ±1.02 


10 




24 ±2 


0.66 ±0.33 


2.50 


5.9 ±3.0 


1.01 ±0.50 


11 




23 ±2 


2.6±1.3 


2.56 


23.4± 11.7 


5.61 ±2.80 


12 




21 ±2 


0.93 ±0.47 


2.75 


8.4 ±4.2 


1.40 ±0.70 


13 




21 ±2 


1.1 ±0.6 


2.81 


9.9 ±5.0 


2.43 ±1.21 


14 




25 ±2 


1.3 ±0.6 


2.03 


11.7±5.9 


2.31 ±1.15 


15 




23 ±2 


0.98 ±0.49 


2.19 


8.8 ±4.4 


2.09 ±1.04 


16 c 




24 ±2 


0.23 ±0.12 


2.50 


20.7 ± 10.4 


0.36 ±0.18 




Aperture radii = 


114" = 6 x (160 [im PSF FWHM) £ 260 pc 


3 




23 ±2 


4.7±2.3 


2.56 


42.3 ±21.2 


11.7±5.9 


7 




24 ±2 


4.7±2.3 


2.28 


42.3 ±21.2 


11.1 ±5.5 


8 




25 ±2 


4.7±2.3 


1.91 


42.3 ±21.2 


10.2 ±5.1 


11 




23±2 


6.7±3.3 


2.62 


60.3 ±30.2 


13.6 ±6.8 


14 




25 ±2 


2.8±1.4 


2.25 


25.2 ±12.6 


4.91 ±2.4 


Total Galaxy 


NGC 6822 


23±2 


83 ± 42 


2.56 


750 ± 380 


103 ±52 



a See Table 1 for source coordinates and cross-identifications. 

k Average dust temperature and implied dust mass over the aperture, derived 
by blackbody fitting. 

c a and dust m ass values derived using the semiempirical SED models of 
Dale ct al. (2001) and Dale & Hclou 12002); typical errors on the a index are 
~ 0.15 (see discussion in 

^Dust mass deri ved using the SED models of lDraine 1HH20ull) . ILi fc Dram! 
I200U|2001) . and lDraine fc Lili2006l) . 

e Note that the local background surrounding region #16 in the IR images is 
very uncertain; see § 13.1.11 



between H i, Ha and IR emission are especially strong 
(though relative variations are evident). 

To further quantify the correlation between H I and 
FIR emission, Figures [7| and [S] show pixel-by-pixel his- 
tograms and scatter plot comparisons of H I column den- 
sity vs. surface brightness measurements in the Ha and 
IR bands. These plots were created by comparing H I, 
Ha, and IR surface brightnesses on an individual pixel 
basis. They clearly demonstrate that FIR and Ha emis- 
sion are strongly peaked in regions with H I column den- 
sities J> 10 21 cm -2 . The 70 and 160 /im histograms show 
that the regions of strongest dust emission (as traced by 
very high-surface brightness 70 and 160 /^m emission; 
thin gray lines in Figure |7J| are associated with the high- 
est H I column densities: the histogram at 16 MJysr -1 
shows a median value ~ 50% higher than the median at 4 
MJysr . There is little difference in the pixel- by-pixel 
distribution of surface brightnesses at Ha or in the IR 
bands; above H I columns of ^5 x 10 20 , the nebular and 
dust emission both increase in strength as the H I col- 
umn increases, with the highest-surface brightness emis- 
sion strongly peaked in regions with H I column densities 
^10 21 cm" 2 . 

We now consider the origin of the evident correlation 
between nebular and IR surface brightness as a function 



of H I column density. If the strength of the dust emis- 
sion is correlated primarily with dust heating, then the 
far-IR colors should become warmer as the ratio of dust 
to H I emission increases. However, if the strength of the 
dust emission is related to an increase in the density of 
dust, then the ratio of infrared to H I emission will not 
be correlated with the far-infrared colors. To investigate 
this scenario, we plot in Figure EH the ratios of S70/SHI 
and Si6o/Shi for the regions listed in Tabled eacn as a 
function of IR color [here, Log(S7o/Si6o)]- These plots 
allow us to gauge how much IR emission is arising per 
unit H 1 mass at a given dust temperature (or, similarly, 
at a given radiation field strength level). It is clear from 
Figure^a) that as the IR color increases (radiation field 
strength rising), the amount of 70 ^m emission per unit 
H 1 mass also increases. At 160 /im, this correlation is 
less evident, suggesting that 70 ^m emission is more sen- 
sitive to local heating, as one would expect for thermal 
emission in the tempera ture range of interstellar dust 
(see also lDale et al.ll2005^ . This arises in part due to the 
70 ^m band falling near the transition point between the 
hot dust emission that dominates shorter wavelengths 
and the ~20-30 K dust emission that radiates at longer 
wavelengths. While the errorbars in both plots are large, 
we can infer that localized heating of dust, in addition to 
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i r 




10 100 10 100 

Wavelength (/zm) 

Fig. 10. — Modified blackbody and SED model fits to the observed data for apertures Al— A8 (see Table 111 and Figure Each panel 
shows the observed infrared flux densities, overlaid with: a single-temperature modified blackbody fit to the observed 70 and 160 fira flux 
densities (fits shown in black); the SED models of Dale & Helou (2002), fitting to the observed 24, 70 and 160 (im flux densities (fits 
shown in blue); the models of Li &c Draine (2001, 2002) and Draine & Li (2006), fitting to all observed flux densities longward of 5 /im 
(fits shown in red); and a stellar component based on population synthesis models, fit to the 3.6 and 4.5 (im bands (green line). Note that 
single-temperature modified blackbody fits strongly underestimate the observed 24 ^m flux densities in all cases. The value of a from the 
Dale & Helou (2002) fits is shown at the bottom right of each panel. 



increased H I column density, plays an important role in 
controlling the nature of IR emission. The different dy- 
namic ranges of nebular, dust and H I emission may also 
suggest increasing H I columns due to photodissociation 
of H2 in regions of active star formation. 

The simplest physical scenario that explains these ob- 
servational trends is one where the H 1 and dust are 
mixed in the ISM, and that as the local radiation field 
strength increas es, so does the to tal IR emission per unit 
dust mass (e.g.. iDale et al.l ll9991. If the increase in the 
dust emission surface brightness was due solely to an 
increase in the local H 1 column density (i.e., a rela- 
tively constant-strength interstellar radiation field pro- 
ducing higher surface brightnesses in regions of larger H 1 
columns) , then the dynamic range of the dust emission 
would be expected to be similar to that in H 1. How- 



ever, Figures and |H1 clearly show a large scatter of IR 
flux densities at a given H 1 column or a given radia- 
tion field strength. We caution that with the present 
observations, it is not possible to assure that the H I, 
Ha and IR observations are all sampling interstellar me- 
dia that are exactly co-spatial; however, the clear cor- 
relations between these various phases (see Figures 0J [7] 
and|SJ suggest that dust absorption and re-radiation are 
closely tied to local star formation intensity. Note that no 
regions with high IR surface brightness and high H 1 col- 
umn density are Ha-deficient (i.e., non-detections); this 
again supports the theory that local star formation is an 
important mechanism in the production of the observed 
IR SED. 

3.2.2. Derivation of Dust Masses 
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i r 




10 100 10 100 

Wavelength (/zm) 

Fig. 11. — Same as Figure UHl but for apertures A9— A16. 



In the following, we calculate the dust mass of 
NGC 6822 (both globally and within the apertures shown 
in Figure |3J) using three different methods of increas- 
ing sophistication. We begin with a single-temperature 
fit, and discuss explicitly the shortcomings of such an 
approach in accurately representing the observed SEDs 
within NGC 6822. We then discuss two SED modeling 
techniques that offer a more physical representation of 
the data. 

First, we apply a simple single-tem perature mo dified 
blackbody fit to th e data (see, e.g.. iHildebrand 119831 
Bia nchi et all 11999ft . Dust mass values are derived via 
the following equation: 

D 2 f 

Mdust = ^B{T) (3) 

where D is the distance to the galaxy, /„ is the flux 
density, k„ represents the absorption opacity of the dust 
(given in Li & Draine 2001), and B(T) is the Planck 
function evaluated at temperature T, derived by fitting 



a blackbody modified by a A 2 emissivity function to 
the 70 and 160 fim data. Using this method, we derive a 
total dust mass of (8.3 ± 4.2) x 10 4 M©; Table H shows 
the inferred dust masses for each region. 

This simple prescription assumes that all the IR emis- 
sion in a galaxy is at both the same temperature and 
in thermal equilibrium. A combination of the composite 
nature of IR emission [arising from large grains, poly- 
cyclic aromatic hydrocarbons (PAH), very small grains, 
etc.], the stochastic heating processes that are of differ- 
ing importance for various types of grains (see Draine 
& Li 2001 for a detailed discussion), and the complex 
relations between dust emission, gas density, and the lo- 
cal interstellar radiation field,, suggests that this single- 
temperature assumption is likely incorrect, at least on 
local scales. Indeed, Figures ^JJ ^2 and ^| show that 
the simple blackbody fits to the 70 and 160 /jm data 
are not intended to represent the hot dust emission at 
shorter wavelengths; the fits fall short of the observed 
24 /im data points by orders of magnitude. Given these 
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Fig. 12. — Same as Figurc^l but for the apertures with expanded radii (see Table|2j and for the global SED of NGC 6822. Note that 
the y-axis has been increased compared to Figures llUI and ITTI 



failures, we seek a more detailed treatment of the SEDs 
that accounts for the potential contribution of multiple 
grain populations to the dust mass. 

A second and more robust method with which to cal- 
culate the dust mass using IR to su b-millimeter da ta is 
with the SEP models presented bv lDale et al.1 i|20flll) and 
iDale fc Heloul <)2002ft . These models assume that large 
grains, very small grains, and PAH molecules contribute 
to the integrated IR energy budget in a galaxy. Semiem- 
pirical SEDs are constructed assuming a power-law dis- 
tribution of dust mass as a function of the strength of the 
local interstellar radiation field (U). With U normalized 
to unity for the radiation field near the sun, masses are 
computed via the relation: 

dM dust (U)ocU- a dU (4) 

where Mdust is the dust mass heated by a radiation 
field of strength < U, and a represents the relative con- 
tributions of the different local SEDs. The fitted val- 
ues of a correlate with an empirical measure of the 



ratio Log\f v (6Qp,m)/.f 1J ( ipQp,Tn)] , which was calibrated 
again st IRAS data (see IDale etalj 120011 IDale k Heloul 
2002). This ratio then corresponds to a scaling factor by 
which the dust mass estimates using single-temperature 
blackbody fits (see above) should be increased to account 
for mass components not well-constrained by the black- 
body fit (see application of this method to NGC 7331 in 
Regan etal. 2004). Values of a between ~ 2-4 imply 
a mass increase of a factor of 9 (see Figure 6 of Dale 
& Helou 2002); a values less than this show a smaller 
increase in dust mass (e.g., a — 1.75 corresponds to a 
factor of ~ 8 increase). 

Fitting these models to the flux densities of the regions 
and the entire galaxy (see Tables and yields values of 
2.0 < a <;2.8 for nearly all regions; Monte Carlo simula- 
tions of the errors on the index a were perfor med during 
the fit ting of these models to 68 galaxies in ID ale et all 
( 2005), providing typical errors of ~ 0.15. The observed 
values of a in NGC 6822 are consistent with those seen 
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in the Dale et al. sample; the only exception is region 8 
(Hubble V), the strongest IR, Ha and radio continuum 
source in the galaxy. This active star formation region 
may have a strong enough local radiation field that a 
large component of cool grains (which would lead to a 
pronounced mass increase) is absent. Figures ITUl 1111 and 
1 121 show these model fits to each region, and Table0]gives 
the values of a for each fit and the inferred dust mass 
using these SED models. Given the uncertainties in the 
dust masses derived from the blackbody fits (see above 
and TableQJ and the varying quality of the least-squares 
fit for each region (compare, e.g., the fits for regions 6 
and 14), all regions are consistent with an increase in 
mass by a factor of ~ 9; the global dust mass increases 
accordingly to (7.5 ± 3.8) x 10 5 M . 

A third and yet more sophisticated estimate of the 
dust mass in the galaxy can be fo und using the mod - 
els of Li & Draine (2001, 2002) and Ipraine fe Lil f20061. 
Here, radiation field strengths are varied via power-law 
distributions; PAH, silicate and graphite grains are il- 
luminated and the resulting SEDs can be compared to 
the observations. The fit of these models to our obser- 
vations is also shown in Figures I10H12I The observed 
SEDs are well-reproduced by these models, when a small 
fraction (~5-30%, depending on the individual region) 
of the dust mass is exposed to radiation fields largely 
in excess of the local average (see Draine & Li 2006). 
The implied dust masses are given in Tables 01 and ^ 
these mas s es are sm aller than those deriv ed using the 
IDale et all l|200lD and lDale fc Heloul l|2002ft models. The 
derived global dust mass is (1.03 ± 0.52) x 10 5 M Q . 

We note that our global flux densities produce SEDs 
that are in good a greement with the fluxes measured in 
llsrael et all ()1996| ) Our use o f soph istica ted SED mod- 
eling techniques [Dale et aD l)2001|) and IDale fe Heloul 
( 2002); Li & Draine (2001, 2002) and Draine & Li (2006)] 
results i n dust mas s es tha t are larger than those de- 
rived bv llsrael et alJ 1^)96). These differences are thus 
attributable solely to our different (and likely more phys- 
ical) modeling approach. 

Note that our data may remain insensitive to a poten- 
tial component of very cold dust (which may be located 
away from regions of star formation). Previous studies 
(e.g., Eales etal. 1989, Devereux & Young 1990, Kwan 
etal. 1992, and Calzetti etal. 2000) have shown that 
observations from the mid-IR into the sub-millimeter 
regime are needed to fully quantify the range of dust 
temperatures, grain sizes, and opacities, and to provide 
a robust estimate o f the total dust mass in a galaxy. 
iPopescu et all l|2002D . for example, noted that ~ 15-20 K 
dust may be present in many nearby galaxies. Data long- 
ward of 160 fim would be needed to determine whether 
such dust is present in NGC6822. Some observational 
results have also demonstrated that dust with tempera- 
tures as cold as 5-10 K (e.g., Madden et al. 2002; Galliano 
et al. 2003, 2005) may be present in dwarf galaxies. How- 
ever, theoretical r esearch (e.g., Li 2004) and other ob- 
servational results Ipumke et alJl2004llBendo et aLll2006L 
e.g.) have suggested that 5-10 K dust temperatures 
would require physically implausible dust grain proper- 
ties or unrealistically high dust masses. Nonetheless, if 
dust cooler than 20 K exists in NGC6822, then these 
data may afford only a lower limit on the total dust con- 
tent of the system. 



3.2.3. The Empirical Dust-to-H I Ratio 

The Spitzer images of NGC 6822 can be compared to 
our sensitive H I imaging to study the dust-to-gas ra- 
tio at ~ 130 pc resolution. We define the dust-to-gas 
ratio (hereafter represented as D) as the ratio of dust 
mass (as traced by these Spitzer observations) to H I 
mass; no correction is applied for the (uncertain) con- 
tribution of CO or H2 gas. At sub-solar metallicities 
(recall that NGC 6822 has a metallicity ~30% Z Q ), the 
molecular phase of the ISM is different from that in 
more metal-rich systems. Theoretical pre dictions (e.g. 
iMalonev fc BlacUl988UBolatto et al.lll,999l) and observa- 
tional evidence re.g..lTavlor et al.l1 998tlLerov etall2005fl 
have shown that CO clouds have a smaller volume fill- 
ing fa c tor in metal-poor systems than in spiral galaxies. 
Ilsraell (|1997t) uses CO observations of selected regions in 
NGC 6822 to infer a total molecular mass ~10% of the 
H 1 mass (though these numbers are highly uncertain, 
depending on the conversion factor used to infer total 
H2 content from observed CO intensities). Since a com- 
plete census of the molecular material in NGC 6822 is 
difficult, we simply study the relation of dust and H 1 
gas; any molecular material in NGC 6822 will lower the 
derived T>. 

Throughout the disk, the mean value of T> (in apertures 
of physical radius ^130 pc, using the dust masses derived 
from the models of Li & Draine) is ^0.004 (with little 
change when using apertures with 4 times larger area; 
see Tables [3] and 0J) . Variations of a factor ~3 are seen 
among individual regions, and the errors on the T> values 
are also high (arising from the relatively uncertain dust 
mass estimates). These average values are lower than 
the average global values found for other SINGS galaxies 
using similar modeling techniques. For example, Draine 
etal. (2006, in preparation) present detailed dust mass 
estimates for many SINGS galaxies: for 40 systems with 
Spitzer data but no submillimcter observations, the av- 
erage global T> (excluding dwarf galaxies and the contri- 
bution of molecular hydrogen) is Mdust/M hi=0.02; i.e. 
a factor of ^5 higher than the average local ratios found 
u sing the same mod els in NGC 6822. 

Israel etal. (1996) finds a global V of 1.4 x 10~ 4 in 
NGC 6822 using IRAS data; our derived global value (us- 
ing the Li & Draine models) is ^8 x 10~ 4 (see Table 
Our Spitzer dat a recover very simi lar global flux densities 
as found in the llsrael et al.l l)1996|) study; the differences 
in dust mass and T> are attributed to the use of detailed 
SED models of the IR emission. Note that all of the 
V values derived for the apertures listed in Table 0] are 
higher than the global value. This is expected since the 
dust is strongly peaked in regions of high H 1 column and 
star formation (e.g., see Figures 03 and [SJ, but the bulk 
of the extended H 1 distribution has little associated dust 
emission (see Figure 0). Recall, however, that our data 
remain insensitive to a potential component of very cold 
dust (see more detailed discussion above). 

We note by examining the individual panels of Fig- 
ures UHl EH and HI that the model fits to the IR SEDs 
predict significant power in the PAH bands between ^5 
and 20 /im (though sec discussion of potential metallicity 
effects in Dale etal. 2005). Most regions show a rise in 
the SED moving from 4.5 /im to 5.8 and 8.0 /im, part 
or all of which may be associated with aromatic band 
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emission in the ISM. Given the metallicity of NGC 6822 
(Z ~ 30% Zq) and th e empirical dependence of PAH 
emission on metallicity En gelbracht et all l)2005|) . it will 
be fruitful to compare the IRAC images with IRS spec- 
tra of various regions in NGC 6822. Initial results of 
SINGS extra-nuclear IRS spectra of the major star for- 
mation complexes (e.g., Hubble I /III, V, X) show that 
PAH emission is prominent throughout the galaxy. This 
will allow a proper spatial and spectral decomposition of 
PAH emission in the metal-poor ISM at unprecedented 
resolution. 

3.3. The Radio - Far-IR Correlation 

In normal star-forming galaxies, there is a remarkably 
tight correlation between the total IR luminosity and the 
strength of radio continuum emission (see, e.g., de Jong 
etal. 1985, Helou etal. 1985, Condon 1992, and refer- 
ences therein). Given the differing physical mechanisms 
giving rise to far-IR and radio emission, the tightness of 
the relation (scatter of less than 0.3 dex over multiple 
order s of m agnitude in luminosity) is especially remark- 
able. iBelll l)2003f) shows that the relation holds for dwarf 
galaxies, which are generally deficient both in IR and ra- 
dio emission with respect to massive galaxies, when the 
emission is normalized by the star formation rate. How- 
ever, these systems require a complex interplay of star 
formation, dust content and magnetic fields to produce 
the apparently linear relation in this regime. 

Figure|3Ih) shows that only sources 2, 3, 4, 7, 8, and 14 
are detected at high S/N in our radio continuum imag- 
ing. These are the six highest-surface brightness Ha re- 
gions in the galaxy. We tentatively detect lower surface 
brightness regions throughout the disk, but the presence 
of strong background sources precludes us from analyz- 
ing them with confidence. Flux densities for these re- 
gions are quoted as upper limits in Table [21 Assum- 
ing that the monochromatic or TIR fluxes and the radio 
continuum emission are related via constant values [q oc 
log(Si R /S RC ); see Bell (2003) and Murphy etal. (2006a) 
for details] , we calculate two values of "q" . The first is a 
monochromatic comparison of the IR and radio contin- 
uum flux densities: 
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The second compares a measure of the total infrared flux 
with the radio continuum: 
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We calculate "q" (see Tabled in individual regions, as 
well as throughout the total galaxy. The resulting values 
are in excellent agreement with those found both locally 
and global ly in larger samples of galaxies. For example, 
IBelll (|2003fl finds a median global qttr = 2.64 ± 02 in a 
sample of more than 200 galaxies. iMurphv etal] l)2006aj) 
explores the radio-IR correlation on a spatially resolved 
basis (sampling physical regions of ^300-750 pc) using 
SINGS data on four nearby spiral galaxies, finding ranges 
of 0.85 <q 24 < 1.10 and 1.95 < q 70 < 2.25 (where q 24 and 
q70 are the monochromatic "q" parameters derived using 
equation [SJ ■ In NGC 6822, there is slight evidence that 
regions with a high Ha/TIR ratio show lower "q" values; 



there is no evidence for a trend of "q" with T> or with 
TIR flux. 

The agreement of the local and global values with the 
canonical radio-IR relations derived for spiral galaxies 
suggests that both the radio and IR are significantly de- 
pleted in NGC 6822, conspiring to produce the observed 
"q" parameters. Different physical mechanisms appear 
to be at work in the ISM of this dwarf galaxy compared 
to those operating in the ISM of more massive spirals. 
In NGC 6822, the bulk of the radio emission appears 
to be thermal: in § 13.1.11 we showed that most regions 
have thermal fractions >90% (the strongest nonthermal 
sources are regions 7 and 8, with thermal fractions of 
^50% and 70%, respectively). In contrast, nonthermal 
processes produce the bulk of radio continuum emission 
in massive spiral galaxies [see Murphy et al. (2006a) and 
Murphy et al. (2006b) for more detailed discussions of 
the radio-FIR correlation in SINGS galaxies] . Evidently 
in NGC 6822 the IR emission per unit star formation is 
reduced by about the ratio of free-free/nonthermal emis- 
sion in normal spirals, conspiring to produce the observed 
"q" values. 

The relative ratios of radio and IR emission strengths 
in NGC 6822 are in marked contrast to the strong de- 
viations from the canonical radio-IR relation found for 
very active low-m e tallicity gal axies in the SINGS sam- 
ple. iCannon et all l|2005l 12006ft show that the radio-FIR 
correlation breaks down in very strongly star forming 
regions. The "supergiant shell" region of IC 2574 has in- 
duced strong variations of the "q" parameter (over an or- 
der of magnitude); the dwarf starburst galaxy NGC 1705 
shows strong FIR emission but is a radio non-detection. 
Other types of galaxies can also be extremely radio- 
deficient, though likely for different reasons than those 
seen in NGC 6822 [see, e.g., the Roussel etal. (2006) 
study of the nascent starburst in the lenticular galaxy 
NGC 1377]. 

These results suggest that star formation intensity has 
a dramatic effect both on the appearance of radio con- 
tinuum emission, and on the derived correlation between 
FIR emission and radio continuum luminosity, in galax- 
ies with small potential wells (for comparison, note that 
more massive starburst galaxies in the Bell 2003 sam- 
ple show little deviation from the derived radio-TIR re- 
lation). The intensity of star formation relative to the 
presence of a stable large-scale disk is important, since 
the ratio of activity to disk density and gravity will gov- 
ern the amount of disruption and blow-out that takes 
place (thus, the irregular nature of NGC 6822 is likely 
also important in determining the ratio of IR to radio 
emission). The relatively uniform values of "q" found 
in larger spiral galaxies and quiescent dwarfs argues that 
disk disruption by star formation is minimal; the strongly 
varying values of "q" in starbursting dwarfs suggests that 
the small disks are being overwhelmed by the star for- 
mation activity. 

4. DISCUSSION AND CONCLUSIONS 

We have presented IR observations of the Local Group 
dwarf irregular galaxy NGC 6822 obtained with Spitzer 
as part of SINGS. The galaxy is highly resolved at all 
imaging bands; the resolution limit is driven by the 
FWHM of the 160 ^m PSF. At ~130 pc physical resolu- 
tion, we study the nature of FIR emission and compare 
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to observations in the optical and radio. These sensitive 
data reveal a wealth of structure within the galaxy. Our 
study confirms some previous results, while also offering 
new insights into the nature of infrared emission in the 
metal-poor ISM. We discuss these topics in more detail 
below. 

The total monochromatic FIR flux densities are found 
to be S 24 = 2.51 ±0.50 Jy; S 70 = 53.2 ±15 Jy; Si 60 = 
136.2 ±40 Jy. SED model fits to these data allow com- 
parison with previous IRAS flux measurements; our 
Spitzer flux densitie s are in ex c ellent agreement with 
the measurements in [Igr^gl^^d.1 1^961) . Using the rela- 
tions presented in[Dalc & Hclou (2002J), these global flux 
densities correspond to a total IR flux of ~ 5.7 x 10~ 12 
Wm" 2 (~ 4.3 x 10 7 Lq). 

IR emission at 24, 70 and 160 /xm is found only within 
the central, high-H I column density (<; 10 21 cm~ 2 ) re- 
gion of NGC6822. The total H I distribution is much 
more extended than the high-surface brightness stellar 
or dust emission components. The 24 /jm surface bright- 
ness profile closely traces that of Ha, delineating regions 
of ongoing star formation. The profiles of 70 and 160 /im 
emission also trace local H I column density maxima. 

Comparing (monochromatic and total) FIR emission 
with a sensitive Ha image reveals a strong morphological 
correlation between IR and Ha surface brightness; here 
we are directly observing the re-processing of photons 
by dust. Roughly 50% of the TIR flux is attributable 
to discrete FIR sources, all of which have at least a 
low- level Ha counterpart; the remaining FIR emission 
arises from a "diffuse" component of the ISM (slightly 
stronger at 70 and 160 fim than at 24 /im; see detailed 
discussion in § I3.1.1J) that is apparently not assoc iated 
with massive star formation. [Galla gher et al.l i|1991f) and 
llsrael et alJ (|1996f) find similar results from IRAS imag- 
ing of NGC6822; we interpret this as the re-processing 
of non-ionizing photons in the ISM or as the escape of 
radiation from the star formation regions. 

We find variations in the relative ratios of nebular 
and dust continuum emission. While differential extinc- 
tion can cause some of the observed fluctuations, it can- 
not account for the full breadth of properties seen in 
the galaxy. Using SINGS observations of the starburst- 
ing dwarf ga laxies IC 2574 and NGC 1705. lCannon etalJ 
l)2005l l2006|) found similar variations and attributed 
them to the direct impact of recent star formation on 
the local ISM (i.e, production and re-processing of pho- 
tons). The results presented here extend this trend to 
more typical, quiescent dwarfs. 

Evidence is mounting from spatially resolved Spitzer 
observations of dwarf galaxies that these systems follow 
a different trend of 24 /Ltm vs. Ha luminosity than what 
is seen in more metal- rich systems (see, e.g., Calzetti 
etal. 2005). This result is in agree ment with previou s 
observations of dwarf galaxies (e.g.. IHunter etall ll989). 
Two different scenarios can account for this observational 
trend. The first one, physically simple and more com- 
monly used, is related to the low dust content of low- 
metallicity galaxies. Here, low-metallicity galaxies sim- 
ply do not have enough metals to form dust in the same 
fraction as more metal rich environments. The second 
scenario posits that dwarfs harbor massive components 
of cold dust (T ~ 5-10 K) that is simply not heated by 
ongoing star formation. Such dust would peak past 200 



/xm and remain undetected in the present Spitzer obser- 
vations, contributing to part of the observed discrepancy; 
observational evidence exists both for and against such a 
component in nearby galaxies (see detailed discussion in 
S I3.2.2JI . Clearly, sensitive observations of a large sample 
of dwarf galaxies between 200 and 1000 /im are needed 
to address this issue. 

As mentioned previously, our global flux densities 
for NGC 6822 are in good agreement with those from 
llsrael et alJ lll996l). We use s ophis ticated SED modelin g 
techniques [Dale etal] (|2001^ and lDale k. Heloul l|2002|) : 
Li & Draine (2001, 2002) and Draine & Li (2006)] to 
de rive dust m a sses t hat are higher than those derived 
bv llsrael et al.l l)1996f) . The differences in mass are thus 
attributable solely to our more physical modeling ap- 
proach, which takes into account the variety of dust 
species and heating levels in each region within the disk 
(see detailed discussion in § I3.2.2fl . For a given TIR lu- 
minosity, our models recover a higher dust mass than 
simple blackbody fits (see detailed discussion in § I3.2.2|l . 
The treatment of multiple grain populations likely means 
that the increased dust masses are a better representa- 
tion of the mass radiating between 10 and 160 /im. Using 
these techniques, we derive a total dust mass of — 10 5 M© 
(this value is derived using the Li & Draine models). 

Comparing to our high-resolution H I imaging, we can 
explore the mass ratio of warm dust to H I gas (2?) in 
the ISM. Given the low nebular abundance of the sys- 
tem (Z ~ 30% Z Q ), we can thus provide constraints 
on a potential trend of T> with abundance. Metal-poor 
dwarf galaxies might be expected to have a lower T> (i.e., 
higher gas mass fraction) th an more metal-rich galaxie s 
such as the Milky Way (e.g ., iThronson & Telescollfl'S^ . 
iLisenfeld fc Ferraral l)1998fl have argued that dwarf ir- 
regular galaxies have a clear relation between these two 
quantities in a global sense (using dust masses derived 
from single-temperature blackbody fits). 

Here, we derive a global V ~ 8 x 10~ 4 . This is -25 
times lower than the average global values found in 
the new SINGS study of Draine etal. (2006, in prep.), 
wherein dust masses are calculated for a large sample of 
galaxies (excluding dwarfs). On local scales (~130 pc), 
the dust mass per unit H I mass is higher (by a factor of 
~5), though still lower than the global values found by 
Draine et al. 

We propose a simple model for the relative ratios of 
dust mass to H I mass in the ISM. As noted above, dwarf 
galaxies may contain large reservoirs of cold dust (T — 
5-10 K) that are undetectable with Spi tzer observations 
(e.g. !MaddeJ27)0llGlmTano et al.ll2005h . Note, however, 
that the dilution of the radiation field required for this 
dust to remain at such low temperatures may require it to 
be located at large distances from the star formation re- 
gions and thus essentially decoupled from the active ISM 
(note that other mechanisms may also be important in 
dwarfs, including dust destruction by supernova shocks; 
see, e.g., Bot etal. 2004). Superposed on this putative 
cool component are regions radiating between 10 and 160 
/im; here, photons are being absorbed and re-radiated as 
IR photons. The low dust contents of metal-poor galax- 
ies plays an important role in the microphysics of this 
process. Galaxies with widespread, elevated star forma- 
tion will heat more of their total dust component than 
more quiescent systems. When more sensitive observa- 
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tions are possible between 200 and 1000 /im, we will be 
able to quantify the relative masses of these components 
in low-mass dwarf galaxies. 

Finally, it is interesting that NGC 6822 seems to fall 
near the canonical measures of the ra dio-FIR relation 
(e.g., |M3 1207)3 iMurohv et all l2006at) . 20cm radio 
continuum imaging reveals six high surface brightness 
sources, each of which appears to be of thermal origin. 
The monochromatic and total IR "q" parameters are 
similar to those found for more massive galaxies. This 
agreement is interesting, given that nonthermal emission 
dominates the radio continua of spiral disks. The sim- 
plest explanation is that depleted radio and far-IR emis- 
sion strengths conspire to produce the observed values. 
These findings in NGC 6822 are in marked contrast to re- 
cent results for very active star-forming dwarf galaxies, 
where significant d eviations from the canonica l radio-FIR 
relation are found ijCannon et allEM 12006]) . Taken to- 



gether, these data suggest that star formation intensity 
is one of the important parameters that governs the ob- 
served radio-FIR relation in dwarf galaxies. 
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